Lattice structures of Larkin-Ovchinnikov-Fulde - Ferrell (LOFF) state 
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Starting from the Ginzburg-Landau free energy describing the normal state to Larkin- 
Ovchinnikov-Fulde-Ferrell (LOFF) state transition, we evaluate the free energy of seven most com- 
mon lattice structures such as stripe, square, triangular. Simple Cubic (SC), Face centered Cubic 
(FCC), Body centered Cubic (BCC) and Quasi-crystal (QC). We find that the stripe phase which 
is the original LO state, is the most stable phase. This result maybe relevant to the detection of 
LOFF state in some heavy fermion compounds and the pairing lattice structure of fermions with 
unequal populations in the ECS side of Feshbach resonance in ultra-cold atoms. 



INTRODUCTION 



It is well known that at sufficiently low temperature, 
an electron with spin up is paired with its partner with 
spin down across the Fermi surface to form a Cooper pair 
with total momentum zero and becomes superconduc- 
tor and exhibits superfluid property. This phenomenon 
is well described by Bardeen-Copper-Schrieffer ( BCS ) 
theory. The most favorable condition for paring is when 
spin up and spin down electrons have the same density. 
Now imagine one apply a magnetic field to split the spin 
up and spin down electrons by Zeeman effect and look 
at the response of a superconductor to the Zeeman split- 
ting. For s-wave superconductor, if the Zeeman splitting 
5^, — /i| — /X| is very small compared to the gap, then the 
superconducting state is stable, if it is much larger than 
the gap, the superconducting state will turn into a nor- 
mal state. When Sjiis comparable to the energy gap Aq 
at zero magnetic field, it may becomes non-trivial. It was 
argued by Fulde and Ferrell^, Larkin and Ovchinnikov^ 
about 40 years ago that an in-homogeneous supercon- 
ductor with pairing order parameter oscillating in space 
may be the ground state at a narrow window of Zee- 
man splitting Sjjii - A0/V2 < < (5/^2 0.754Ar,3-'^ ( 
Fig.l ). This in-homogeneous state is called LOFF state 
where the Cooper pairs carry a finite momentum. In FF 
state, A(a;) = Ape*^'^ where q ^ kp^ — kp\^, the Cooper 
pairs carry finite superfluid momentum , while in the LO 
state, A(a;) = Aocosq* • x, the Cooper pairs carry two 
opposite momenta. The LOFF state breaks both C/(l) 
gauge symmetry and translational order. Unfortunately, 
so far, the LOFF state has never been observed in con- 
ventional superconductors, because in these systems, the 
Zeeman effect is overwhelmed by orbital effects. How- 
ever, this LOFF state has attracted renewed interests 
in the context of organic, heavy fermion and high 
cupratea^, because these new classes of superconductors 
may provide favorable conditions to realize the LOFF 
state. Recently, experiments^ on penetration depth mea- 
surement on CoC'eln^ shows that at a temperature be- 
low 250 mK, for magnetic field applied parallel to the ab 
plane, two phase transitions were detected, one of which 
maybe identified as a phase transition from LOFF state 
to normal state transition. Also the measurement of ther- 
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FIG. 1: The phase diagram of LOFF state. Sfi is the Zeeman 
splitting, T is the temperature, Ao is the energy gap at the 
balanced case Sjj. — 0. 



mal conductivity^ on CoCeln^ shows anisotropy in real 
space, which could be interpreted as domain wall forma- 
tion, namely, a stripe phase but possibly with higher har- 
monics. LOFF states also played important roles in high 
density quark matter, astrophysics^ and superconductor- 
ferromagnet heterostructures^. With the development of 
trapped cold atoms system, it was proposed that due to 
absence of orbital effects, ultracold neutral fermion gases 
with unequal populations may realize the LOFF state in 
a tiny window on the BCS side of Feshbach resonance^. 
Recently, it was argued inii that the LO state, in fact, 
may be stable in an appreciate regime in the BCS side of 
the Feshbach resonance. 

Before we discuss the phase diagram Fig.l, we re- 
viewed the basic facts of classical Lifschitz point which 
is closely related to normal state to LOFF state phase 
transition. This connection is not that new, but has not 
been stressed in any literature. The free energy near a 
classical {d,d±) Lifshitz point is^^: 

iJ = i y d'^x[tm^ + K\\ (Vm)2 + K^{Wmf 

+ Liy'^mf] +uj d'^xm^ + ■■■ (1) 



where > and m{x) is a n > 2 component order 
parameter, the dimension d is divided into d± perpen- 
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FIG. 2: (a) Phase diagram of Classical Lifshitz point (CLP). P 
is the Paramagnetic phase, F is the ferromagnetic phase, M is the 
modulated phase. The LP point is at {t,K±^) = (0,0). The dashed 
line is the P-M transition we are studying, (b) Momentum shell of 
width A around 2d roton surface. 

dicular dimension and parallel dimension. Its phase 
diagranii^ is shown in Fig. 2. 

Let me review the phase transition from P to M tran- 
sition along the dashed line shown in Fig. 2. In the P 
phase along the path close to the P-M transition bound- 
ary, t > 0, i^Ty > 0, K± < 0, for simplicity, we can set 
— 0, the propagator -D(fc|| = 0, fcj^) can be written 
as D{kj_) = t + K±_k']_ + Lk\ = A + L{k\ - k^f where 

A = t — , k^ = ^^j^ ■ It is easy to see the minima is lo- 
cated at the " roton " surface fc^ ( Fig. 2b) , in sharp con- 
trast to K± > case where the minimum is at fcj^ = 0. 
This class of problems with minima located at /c^ > 
was first investigated ini^ and has wide applications in 
the context of liquid crystals^^. When A > 0, the system 
is in the paramagnetic ( P ) phase with < m >= 0, while 
when A < 0, it is in a modulated ( M ) phase with the 
mean field structure < m{x) >— J^f^i \e^^''^, li = K- 

The P — M transition happens at A = 0, namely, t — -jj^ 
as shown in Fig. 2. The M phase breaks both the inter- 
nal 0(n) rotational symmetry and the translational sym- 
metry, therefore supports two kinds of Goldstone modes: 
phase mode due to the 0{n) symmetry breaking and the 
lattice phonon mode due to the translational symmetry 
breaking. At the mean field theory, the P-M transition 
is 2nd order. Under fluctuations. For d± = 1, the ro- 
ton surface in Fig. 2b, in fact, turns into two isolated 
points, the transition which describes nematic-Smectic 
A transition in liquid crystal remains 2nd order. How- 
ever, for d± > 2, the transition becomes a fluctuation 
driven 1st order transition as shown by Rcnormalization 
group analysis in^^. Indeed, to some extent, the LOFF 
phase diagram Fig. 2 looks similar to Fig. 1 if we identify 
Zeeman splitting Sfj, as the pressure —K±, normal phase 
as the paramagnetic phase, the superconducting phase 
as the ferromagnetic phase and the LOFF state as the 
modulated phase. 

Of course, the original pairing problem of fermions 
with unequal populations are a fermionic problem. How- 
ever, just like usual normal state to BCS superconductor 



transition, one can integrate out fermions at any finite 
temperature and lead to the following Ginsburg-Landau 
free energy describing the normal state to the LOFF state 
transition^ii^iii: 

+ d[irny^r+^'{yrr]+e\ijf, (2) 

where qo kp^ — kp^. 

Indeed, this action is very similar to the Lifshitz ac- 
tion EqnlT]with K±^ < 0, so similar procedures following 
Eqn[l]can be used. Substituting ■0 = X^g^gg*'^^ where 
G are the shortest reciprocal lattice vectors into the above 
equation and combining terms lead to the GL free energy 
in momentum space: 

/ = X^^^'gI^gP +wX!'^Gi'0G2V'G3V'G4'5Gi+G2+G3+G4 

G ^ G 

-1- W ^ -001 V'Gs V'Gs V'G4 V-Gs V'Ge ^Gi -HG2 +G3-(-G4+G5 +(^^?) 
G 

where r = T — Tc and u, v are functions of the coefficients 
&, c, d, e in Eqnl^] and G. 

If r > 0, the system is in the normal state with 
< '0(6') >= 0, while when r < 0, it is in a modulated 
( M ) phase with the mean field structure < 0(a;) >= 
YJi^i A,e'9'-^, qi = qo. This M phase is the LOFF state. 
The LOFF state breaks both U{1) symmetry and the 
translational symmetry, therefore it supports two kinds 
of Goldstone modes. (1) the Goldstone mode due to the 
U{\) symmetry breaking, but it was "eaten" by the gauge 
field due to Higgs mechanism in electron pairing case in 
condensed matter system, but will stay in the neutral 
atom pairing case in ultra cold atom atomic experiments 
(2) the lattice phonon modes due to the translational 
symmetry breaking, they will survive the gauge field fluc- 
tuations. In this paper, we approach the LOFF state 
from the normal state and try to determine what is the 
lowest lattice structure of the LOFF state. P = 1 cor- 
responds to the FF state, P = 2 corresponds to the LO 
state. It is known that the FF state, being carry finite 
superfluid momentum, is always unstable. The LO state 
has nodes where the excess fermions reside. However, it 
is still not know the LO state is the most favorable lattice 
structure. In this paper, we will study what is the low- 
est lattice structure by considering seven most common 
lattice structures namely the stripe, square, triangular. 
Simple Cubic (SC), Face centered Cubic (FCC), Body 
centered Cubic (BCC) and Quasi-crystal (QC) listed in 
Table I. The stripe case corresponds to the original LO 
state. 

The rest of the paper is organized as follows. In sec- 
tion II, we compute the coefficients of the free energy of 
the LOEF states with different lattice structures. In sec- 
tion III, by comparing the free energy and the transition 
temperature of all the seven lattice structures of LOFF 
state, we find the lowest energy lattice structure remains 
the LO state. In the appendix A, we discuss in detail how 



(a) m=2 (b) m=4 (c) m=6 

FIG. 3: The set of shortest reciprocal lattice vectors Lq for one 
and two dimensional lattices (a) Stripe lattice (b) Square lattice 
(c) Triangular lattice 

to get the geometrical factors in the fourth and sixth or- 
der terms which are used in evaluating the free energy 
of the seven lattices. As a byproduct, we corrected some 
over-counting mistakes in describing liquid to solid tran- 
sition in the textbook In appendix B, we revisit the 
solid to liquid transition by considering both cubic and 
quartic term and show that the BCC lattice remains the 
favorable lattice in the presence of cubic term in a certain 
region. 



2. EFFECTIVE FREE ENERGIES OF THE 
LOFF STATE WITH DIFFERENT LATTICE 
STRUCTURES 

We only look at the subset Lq spanned by all the 
shortest reciprocal lattice vectors G — qo- In the ground 
state, tpc has to be real up to a global phase. From 
the point group symmetry of the lattices, ipc is a con- 
stant when G belongs to Lq- FoUowing^'^, we have scaled 
uq — > nQmT^/'^ so the quadratic term is the same for all 
the lattices. Then Eqnl3]is simplified to the effective free 
energy in different lattices: 

/ = \'^^G + "aV'G + '"ai'G (4) 

where a stands for different lattices. In the following, 
we will calculate the fourth order term Ua and the sixth 
order term Va for different lattices respectively. 

1. The fourth order term Ua- For stripe phase, 
square lattice , triangular lattice, SC and FCC, as shown 
in the appendix A, there are only contributions from 
paired vectors to the quartic term = 3(1 — ^) where 
771 is number of the vectors in the set Lq. Therefore 
u\\ — |m, un = |u,MA = |m, Use = ^u,Ufcc = ^u. The 
set Lg for different lattices are shown in Fig. 3 for one and 
two dimensional lattices and Fig. 4 for three dimensional 
lattices. 

For a BCC lattices, there is an additional vertex con- 
tribution Uy = u coming from the 4 vectors from any of 
the six vertices. So in all, u^cc — Up + Uy — ^u. 

For a quasi-crystal, we have an additional contribution 
from the non-planar diamond»l^ Unpd = ^u, so in all, 
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2. The sixth order term For the stripe phase, 
square lattice, SC and FCC, there are only contributions 



(a) m=6 (b) m=12 ((,) ni=8 (d) m=30 

FIG. 4: The set of shortest reciprocal lattice vectors Lq for three 
dimensional lattices (a) Simple Cubic (b) BCC lattice (c) FCC 
lattice (d) Quasicrystal 




(a) (b) (c) 

FIG. 5: non-paired contributions to sixth order term in BCC lat- 
tice (a) a pair of opposite vectors plus four vectors coming out of 
one vortex, lOv, (b) a non-planar triangle diagram with the com- 
mon edge chosen twice, (c) a triangle diagram, each vector in 
the triangle was chosen twice, ; for the triangle lattice in Fig. 3c, 
this term is f n 



from paired vectors — 5(3m^ — 9m + 8)/to^. So we 
get v\\ = 2^v,va = 6jv,Vsc = ^v,Vfcc = lOw. 

For the triangular lattice, there is an additional con- 
tribution vtri = ft! coming from the closed triangles dia- 
gram ( Fig. 5c ). So we get va = Vp + vtri = 9|w. 

For the BCC, in additional to the paired vector contri- 
butions yP = ^^^j there are also contributions coming 
from the three configurations listed in Fig. 5 which is 
In all, Vbcc = 220/9i;. 

For Quasicrystal, in additional to the paired vector 
contributions — ^-^v, there are also contributions 
coming from the four configurations listed in Fig. 6 which 
is ^v. In aU, Vqc = 497/18w 

The Ua and Va for the seven lattices are listed in the 
following table. 
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TABLE I; u and v for the seven lattices 
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difference between Quasicrystal and BCC 
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FIG. 6: non-paired contributions to sixth order term in Quasicrys- 
tal lattice (a) a pair of opposite vectors plus a non-planar diamond 
structure, ^v. (b) a non-planar triangle diagram with the com- 
mon edge chosen twice, ^v. (c) a triangle diagram, each vector 
in the triangle was chosen twice, j^v- (d) two triangles with no 
common edges, 
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FIG. 7: u is positive. Difference between gqc and gtc 



3. OPTIMAL LATTICE STRUCTURE OF THE 
LOFF STATE 

In the original GL action EqnO u can be negative and 
positive. In case v is also negative, then an eighth order 
is needed. In this paper, we assume v is always positive 
to keep the system stable. In the following, we discuss 
u < and u > cases respectively. 

1. u is positive. It is easy to see that U|| < «□ < Ugc = 
UA < Ufcc < Ubcc and v\\ < vq < Vsc < va < vjcc < 
Vbcc so for any given xp: /||(^) < /□(V') < fscW < 
< ffccW < fbccW- Then /||(^||) < /□(V'n) < 
fsc{ipsc) < /a(V'a) < ffccii'fcc) < fbcci'i^bcc)- However, 
more work is needed to compare Quasicrytal with BCC. 
Minimization of Eqn|4] leads to the order parameter and 
the free energy: 
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(5) 



Defining r = where x is dimensionless and plug- 
ging it into EqnO we get fa = ^ga{x) where ga are 
dimensionless functions and a stands for Quasicrytal and 
BCC. Comparing these two functions, we find that there 
is a shift of order between these lattices as shown in Fig. 7. 

When -0.274^ < r < 0, g^^ < g^cc thus /^^ < hcc- 

2 

However when r < -0.274—, g^c > gtcc thus fqc > hcc- 
In any case, the stripe phase is the lowest free energy 
lattice. 

2. u is negative. Eqn. [5] still hold for m < 0. We can use 
the same method used when u is positive. Defining r ~ 
and plugging it into Eqn[5l we still have the following 
expression = ^ga{x). For seven different lattices, we 

get the same coefficient but different functions ga 
with respect to x. 
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FIG. 8: u is negative, (a) gdx) of seven different lattices, it 
is hard to see the difference between FCC and triangular in 
this scale, (b) The difference between triangular and FCC in 
the expanded scale. 



Comparing ga, g\\, gA, gbcc, g/cc, gsc, gqc shown in 
Fig. 8a, we find that there is a shift of order between 
triangle lattice and FCC lattice shwon in Fig.8(b). The 

transition temperature of FCC is Tfc 



and that of triangular lattice is Ta 
It shows that as the temperature is decreased, the 
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first solid phase between these two is FCC, but 
when the temperature is further decreased below 
the transition temperature of triangular lattice and 

2 

when r < —0.617-^, the triangular lattice has the 
lower energy than FCC, which means that FCC is 
a mestable state after that. In general, we have 

2 

the following relations, when —0.617^ < r < T/cc, 
9\\ < 9o < 9sc < 9fcc < 9A < 9bcc < 9qc thus 
f\\ < /□ < fsc < ffcc < /a < fbcc < fqc- When 

2 

r < -0.617—, .g|| < ga < gsc < 9a < 9 fee < 9bcc < 9qc 
thus/|| < /□ < fsc < /a < ffcc < fbcc < fqc- In any 
case, the stripe phase is always the lowest energy state 
of all the seven lattices. 

In fact, we can get the same result from the critical 
transition temperatures of different lattices. It is known 
that the transition temperature in the above model is 

2 

Tc = Plugging Ua and Va for different lattices, we 

find out that the stripe lattice has the highest transition 
temperature as expected, which means when we decrease 
the temperature, the first solid phase will be the stripe 
phase. 



4. CONCLUSIONS 

In this paper, we study the transition from the normal 
state to the LOFF state from the GL free energy in a 
mean field theory. We consider seven most common lat- 
tices. By comparing the free energy and the transition 
temperature of the seven lattice structures, we find that 
the lowest energy lattice structure of the LOFF state is 
the stripe phase, which is the LO state originally pro- 
posed by Larkin and Ovchinnikov- . Our result shows 
that in heavy fermion system or cold atom system, at a 
sufficiently low temperature, if a LOEF state can be re- 
alized, then its lattice structure will likely to be a ( stripe 
) LO phase which will lead to anisotropy in many phys- 
ical measurable quantities. Although so far, there is no 
direct probe on the structure of the order parameter in 
all these heavy fermion materials, in experiment in^, the 
thermal conductivity measurement was used to probe the 
anisotropy of the order parameter, especially the struc- 
ture of the nodes in the momentum space. The experi- 
ment indeed show the anisotropy of the thermal conduc- 
tivity of CoCeln^ in the possible LOFF state regime in 
Fig.l. Our results suggest that the LOFF state observed 
in the experiment is the original LO state. Of course, the 
order parameter may contain higher Harmonics terms. 
Recently, it was argued in— that the LO state may be 
stable in an appreciable regime in the imbalance versus 
detuning phase diagram in the BCS side of the Feshback 
resonance. It is not known if the GL action still can be 
used to describe the normal to the LO transition at T = 
where r = p — pc where Pc is the critical polarization dif- 
ference, because at T = 0, the residual fermions can not 
be integrated out, especially near the transition point. 
However, we expect the normal to the LOFF state tran- 



sition is still of the Lifshitz type first order transition. 
Well inside the LOFF state, mean field analysis in the 
paper still holds, so the results still apply. 

We thank Kun Yang for helpful discussions and Yong 
Tang for technical support. The Research at KITP was 
supported in part by the NSF under grant No. PHY-05- 
51164. 

APPENDIX 



APPENDIX A: NUMERICAL FACTOR IN BCC 
AND QUASICRYTAL LATTICES 

In this appendix we present in detail the procedures 
to get the numerical factors for the forth and sixth or- 
der term used in the main text. There are many ways 
to draw the direction of arrows in the diagrams in the 
main text. Of course, all the different ways should 
give exactly the same numerical factors. But for some 
choices, special cares are needed to avoid overcounting 
the contributionsi^.. In the main text, we just showed 
the most convenient choice. 

Now there are two methods to get the paired vector 
contribution to the forth and sixth order term. The first 
method A is a constructive method by which we count 
the number of ways the i^iG) can take one by one, this 
way is straightforward and can naturally avoid any possi- 
ble over-countings, but it is a little bit tedious especially 
when the order increases. The second method is by some 
combination trick method B, this way is less straightfor- 
ward, but can be more effective when the order increases. 
The agreement of the final results between the two meth- 
ods can insure the correctness of our results. 

Method A: 

For the forth order term, the first ipc can take m 
choices and then (1) the second V'g takes the same vec- 
tor again and then the next two ipc must take exact 
opposite of that vector, so there is only 1 choice here 
(2) the second ipc takes the opposite of the first vec- 
tor. Then for the third and forth ipc have to be opposite 
and have m choices. This case essentially reduces to the 
quadratic case. (3) the second tpc takes one of the m — 2 
choices which is different than the first vector and its 
opposite. Then the third and forth must be the exact 
opposite of the first and second vector, therefore there 
are only 2 choices here. The total sum of all the choices 
are m[l + m + 2(m — 2)] — 3m{m — 1). After rescaling 
by m^, we get uP{m) — 3(1 — ^)u. 

For the sixth order term, the first can take m 
choices and then (1) the second V'g takes the same vector 
again (a) the third -00 also take the same vector, then 
there is only one choice left for the rest three i/'G- (b) 
the third ij^c also take the opposite vector, then from 
the calculations in the forth order term, then there are 
3to — 3 choices for the rest three ipG- (c) the third ■i/'G 
takes one of the m — 2 choices which is different than the 
first vector and its opposite, then there are 3 choices for 



6 



the rest three ijjG- So, adding (a) + (&) + (c), there are 
1 + (3m — 3) + 3(r7i — 2) = 6to — 8 choices for case 1. 
(2) the second ipc takes the opposite of the first vector. 
Then this case essentially reduces to the forth order case, 
so there are 3m(m — 1) choice. (3) the second tAg takes 
one of the m — 2 choices which is different than the first 
vector and its opposite, (a) the third takes one of 
the first two choices, there there are 2x3 choices (b) 
the third V'G takes the opposite of one of the first two 
choices, there there are 2 x 3(m — 1) choices, (c) the 
third ipQ takes one of the to — 4 choices which is different 
than the first and the second vectors and their oppo- 
sites, then there are 6 choices. So, adding (a) + (&) + (c), 
there are (to — 2)[6 + 6(to — 1) + 6(to — 4)] = 12(to — 2)^ 
choices for case 3. Adding all the (1) + (2) + (3) = 
TO[6TO-8 + 3TO(TO-f) + 12(TO-2)2] = 5to(3to2-9to + 8). 
After rescaling by m? , we get w^J — 5(3 — 9/to + 8/to^)w. 
Method B. 

For the forth order term, There are 2 choices: (1) we 
choose same pair twice or choose two different pairs. If 
we choose same pair twice, first we have ^^ choices of 
paired vector, then we put this pair into 4 location. The 

contribution of this is ^2: x = 3to. (2) We choose 

two different pairs, which ( 2 j ^'^'^ them into 4 

different location, that will be 4!. So this term will give 
3to(to — 2). The sum of the above two contributions gives 
3to(to — 1). Rescaling by m? , we get = 3(1 — l/m)u 
which is the same as that achieved by the method A. 

For the sixth order term: (1) we choose the same pair 
three times and put them into 4 locations, which is ^ x 



= f Oto. (2) we have two pairs with one pair chosen 



twice, that will be 2 x ( ^ 
' 2 



'6^ 
X, ,x 



X 2! 45to(to — 



2). (3) we have three different pairs in 4 locations, which 

m \ 

^ j x6! = 15to(to — 2)(to, — 4). The sum of the above 

three contributions gives 5to,(3to^ — 9m + 8). Rescaling 
by m^, we get v^, = 5(3to^ — 9m + 8)/m'^v which is the 
same as that achieved by the method A. 

Next, we are going to show how to get the nontrivial 
terms for BCC and Quasicrytal. 

1. BCC lattice 

For BCC, we can see from Fig. 4 that for each vertex, 
there are two arrows coming in and two arrows coming 
out which, in spin ice case, is called "two in, two out" 
rule and the sum of four vectors from any vertex must 
be equal to 0. 

a. The forth order term: So we have this nontrivial 
vertex contribution to the forth order term. There are 6 
vertices, and therefore there are 6 sets of these vectors. 
Their contribution to the forth order term after rescaling 

is Uyet = u. 

b. The sixth order term: In addition to the paired 
contributions calculated by the two methods above, 



there are three non-paired contributions listed in Fig. 5. 
(5a) there is a paired vectors plus four vectors coming 
out of any of the 6 vertices. There are six pairs of 
vectors and 6 vertices. If we choose any vertice, there are 
4 pairs of vectors we have exactly one vector already 
chosen inside of the vertice So the contribution of this 

6x4x (j x4! = 8640. And also we have 2 pairs 

having no same vector as that in this vertices.. This 
contribution is 6 x 2 x 6! = 8640. The sum of these two 
terms after rescaling gives lOv. 

(5b) two different triangles having a common edge. 
For each edge, we have exactly one of these choices, so 
there are going to be 12. We want to put one triangles 
into 6 locations. The number of way to do that is 
^6^ 

(5c) by observation the sum of three vectors from a 
closed triangles equals to 0, therefore there is a contri- 
bution coming from one closed triangle with each side 
chosen twice and there are 8 different closed triangles. 

/6\ ^4> 

The number of ways to do that is 8 x 



12 X { "] X 4! = 4320. After rescahng,we have |u 



, X , 
2/ V2 



720. 



After rescaling, we have |w 
Note that two triangles having no common edge con- 
tribution has already been included in the (5b) and 
(5c). After the sum of (5a), (5b), (5c) and rescale, we get 
Vbcc = 220/9V. 

2. Quasicrytal lattice 

a. The forth order term: FoUowing-'^^, in addition to 
the paired vector contribution calculated by the methods 
above, there are also 30 non-planar diamond contribu- 
tion(Fig.6(a)) to the forth order term. After rescaling, we 
have Uupd = fu- 

b. The sixth order term: Following the same procedure 
for the BCC lattice, in addition to the paired vectors con- 
tribution, there are four non-paired contributions listed 
in Fig. 6. (6a) one paired vectors plus any non-planar di- 
amond structure. There are 15 pairs of vectors and 30 
non-planar diamond structure. Their contribution after 

rescaling is 30x4x x 4! 4-30 x (15-4) x 6! = 280800. 

After rescaling, we have ^v. 

(6b) we have non-planar closed triangles with a common 
edge. For each edge, there is exactly one such config- 
uration, so there are 30 of them. Their contribution is 
6^ 

(6c) there is closed triangle with each sides chosen twice. 
In Qccrytal, there are 20 different closed triangles. Their 

contribution is 20 x x = 1800. After rescaling, 
we have j^v 

(6d) a contribution from two different triangles with no 
common edges. For each triangles, there are 12 different 
triangles that haven't been included in previous contribu- 
tions.Their contribution is 20 x 12 x 6!/2 = 86400. After 



30 X ( ^ I X 4! = 10800. After rescaling, we have ^v. 
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rescaling, we have ■^v 

After the sum of all these terms plus the trivial contri- 
bution from paired vectors, we get Vqc = 497/18?;. 



APPENDIX B: LIQUID TO SOLID 
TRANSITION, REVISIT 



The liquid to solid transition was studied ini^ by con- 
sidering only the cubic term. In this appendix, we will 
consider the effects of both the cubic and forth order 
term. For liquid to solid transition, expanding the order 
parameter to the forth order term, we have 



+G3,0 



G 



-I- u > 6^ ,^ 



2+G3+G4,0 ■ 



(Bl) 



Obviously, the difference between liquid to solid transi- 
tion and the normal state to LOFF state transition con- 
sidered in the main test is that there is a cubic term in 
the former, but not in the latter. FoUowingi^, one can 
simplify Eqn IBll to: 



(B2) 



Because the Quantum Hall to insulator transition in 
single layer quantum Hall system—, Excitonic super- 
fluid to Excitonic solid transition in electron-hole bilayer 
system^^ happen in two dimensions, we will first compare 
two dimensional lattices, namely stripe lattice, square 
lattice and triangular lattice. For stripe and square lat- 
tice, it is easy to see the cubic term Wa — 0, because 
there is no closed triangle in all these lattices. Minimiz- 
ing the free energy we have /„ = — . Since uy < 
it is very easy to see that /y < /□. 

For triangular lattice, the contribution to the cubic term 
from a closed triangle was evaluated ini^ to be wa = 
^w. From the appendix A, we get u/\ = 2.5u. Mini- 
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mizing the free energy Eqn lB2l leads to: 



-3w 



71A = 



/a = 



Sua 

1 w^r , 

r — 

4 4ma 4ma 



riA (r I — 

4^ 4ma 



4m 



A 



4u 



A 



(B3) 



We can see although the stripe lattice doesn't have a 
cubic term, but u\\ < ma. The more complete way to 
evaluate which one has a lower free energy must take 
these two terms into consideration, not just considering 
the cubic term as did ini^. Now Define r = —x and 

u 

find out the difference between /A,/y, we compare 
numerically the two functions within the range r < 0. 
We still find that the triangular lattice always has a 



triangular 




FIG. 9: Function / for 3 two dimensional lattices: /□ , /a , /|| 



lower free energy than stripe lattice when it is close to 
the transition point as shown in Fig. 9. When the tem- 
perature is further decreased, we find that numerically 
/a > f\\ ■ But it is known that GL theory is only valid for 
weak first order transition and second order transition 
close to transition point. If the temperature is further 
decreased, the validity of Eqn lB2l mav be questioned. As 
shown in Fig. 9, we still find that the triangular lattice 
always has a lower free energy than square lattice. 

Now we will generalize the above consideration to the 
seven lattices listed in table L For stripe, square, SC and 
FCC, it is easy to see the cubic term Wa = 0, because 
there is no closed triangle in all these lattices. Minimizing 
the free energy we have fa ~ — -^g^ . 

For triangular lattice, BCC and Quasicrystal, the con- 
tribution to the cubic term from a closed triangle was 
evaluated in^^ to be Wa = '^1 \prri after rescaling. Mini- 
mizing the free energy Eqn lB2l leads to: 



_ —Zwa + y^9w^ — IGwar 

= - 

" 4^ 4u„ 4u„Mu„^"" 4^'' AujAua 

(B4) 

Following the same method used previously, we can de- 
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difference between SC and BCC 



difference between Quasicrystai and BCC 




fine r = —x and find the difference between fa- We com- 
pare numerically the functions within the range r < 0. 
Since we know that u^c < w/cc, fsc < ffcc, so in or- 
der to find which one has the lowest energy, we only 
need to compare triangular lattice, SC and Quasicrys- 
tai with BCC. In Fig. 10, we only show the difference 
between BCC, Quasicrystai, triangular and SC. 



(a) 



(b) 



difference between trianguiar and BCC 




FIG. 10: Difference between (a) fsc and ftcc, (b) fqc and fbcc 
(c) /a and ftcc 



We find that when the temperature is decreased just 
below the transition temperature of the lattices, the lat- 
tices with a cubic term have a smaller free energy than 
the lattices which do not. Fig. 10 shows that BCC lat- 
tice has the lowest free energy and the highest transition 
temperature in a range just below the transition point. 
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